Do quasi-regular structures really exist in the solar photosphere? I. Observational evidence

Two series of solar-granulation images -- the La Palma series of 5 June 1993 and the SOHO MDI series of 17--18 January 1997 -- are analysed both qualitatively and quantitatively. New evidence is presented for the existence of long-lived, quasi-regular structures (first reported by Getling and Brandt (2002)), which no longer appear unusual in images averaged over 1--2-h time intervals. Such structures appear as families of light and dark concentric rings or families of light and dark parallel strips (``ridges'' and ``trenches'' in the brightness distributions). In some cases, rings are combined with radial ``spokes'' and can thus form ``web'' patterns. The characteristic width of a ridge or trench is somewhat larger than the typical size of granules. Running-average movies constructed from the series of images are used to seek such structures. An algorithm is developed to obtain, for automatically selected centres, the radial distributions of the azimuthally averaged intensity, which highlight the concentric-ring patterns. We also present a time-averaged granulation image processed with a software package intended for the detection of geological structures in aerospace images. A technique of running-average-based correlations between the brightness variations at various points of the granular field is developed and indications are found for a dynamical link between the emergence and sinking of hot and cool parcels of the solar plasma. In particular, such a correlation analysis confirms our suggestion that granules -- overheated blobs -- may repeatedly emerge on the solar surface. Based on our study, the critical remarks by Rast (2002) on the original paper by Getling and Brandt (2002) can be dismissed.Comment: 21 page, 8 figures; accepted by "Solar Physics

Pathways of Large-scale Magnetic Couplings Between Solar Coronal Events

The high-cadence, comprehensive view of the solar corona by SDO/AIA shows many events that are widely separated in space while occurring close together in time. In some cases, sets of coronal events are evidently causally related, while in many other instances indirect evidence can be found. We present case studies to highlight a variety of coupling processes involved in coronal events. We find that physical linkages between events do occur, but concur with earlier studies that these couplings appear to be crucial to understanding the initiation of major eruptive or explosive phenomena relatively infrequently. We note that the post-eruption reconfiguration timescale of the large-scale corona, estimated from the extreme-ultraviolet afterglow, is on average longer than the mean time between coronal mass ejections (CMEs), so that many CMEs originate from a corona that is still adjusting from a previous event. We argue that the coronal field is intrinsically global: current systems build up over days to months, the relaxation after eruptions continues over many hours, and evolving connections easily span much of a hemisphere. This needs to be reflected in our modeling of the connections from the solar surface into the heliosphere to properly model the solar wind, its perturbations, and the generation and propagation of solar energetic particles. However, the large-scale field cannot be constructed reliably by currently available observational resources. We assess the potential of high-quality observations from beyond Earth's perspective and advanced global modeling to understand the couplings between coronal events in the context of CMEs and solar energetic particle events

Constructing and Characterising Solar Structure Models for Computational Helioseismology

In this paper, we construct background solar models that are stable against convection, by modifying the vertical pressure gradient of Model S (Christensen-Dalsgaard et al., 1996, Science, 272, 1286) relinquishing hydrostatic equilibrium. However, the stabilisation affects the eigenmodes that we wish to remain as close to Model S as possible. In a bid to recover the Model S eigenmodes, we choose to make additional corrections to the sound speed of Model S before stabilisation. No stabilised model can be perfectly solar-like, so we present three stabilised models with slightly different eigenmodes. The models are appropriate to study the f and p1 to p4 modes with spherical harmonic degrees in the range from 400 to 900. Background model CSM has a modified pressure gradient for stabilisation and has eigenfrequencies within 2% of Model S. Model CSM_A has an additional 10% increase in sound speed in the top 1 Mm resulting in eigenfrequencies within 2% of Model S and eigenfunctions that are, in comparison with CSM, closest to those of Model S. Model CSM_B has a 3% decrease in sound speed in the top 5 Mm resulting in eigenfrequencies within 1% of Model S and eigenfunctions that are only marginally adversely affected. These models are useful to study the interaction of solar waves with embedded three-dimensional heterogeneities, such as convective flows and model sunspots. We have also calculated the response of the stabilised models to excitation by random near-surface sources, using simulations of the propagation of linear waves. We find that the simulated power spectra of wave motion are in good agreement with an observed SOHO/MDI power spectrum. Overall, our convectively stabilised background models provide a good basis for quantitative numerical local helioseismology. The models are available for download from http://www.mps.mpg.de/projects/seismo/NA4/.Comment: 35 pages, 23 figures Changed title Updated Figure 1

Acoustic Events in the Solar Atmosphere from Hinode/SOT NFI observations

We investigate the properties of acoustic events (AEs), defined as spatially concentrated and short duration energy flux, in the quiet sun using observations of a 2D field of view (FOV) with high spatial and temporal resolution provided by the Solar Optical Telescope (SOT) onboard \textit{Hinode}. Line profiles of Fe \textsc{i} 557.6 nm were recorded by the Narrow band Filter Imager (NFI) on a $82" \times 82"$ FOV during 75 min with a time step of 28.75 s and 0.08$"$ pixel size. Vertical velocities were computed at three atmospheric levels (80, 130 and 180 km) using the bisector technique allowing the determination of energy flux in the range 3-10 mHz using two complementary methods (Hilbert transform and Fourier power spectra). Horizontal velocities were computed using local correlation tracking (LCT) of continuum intensities providing divergences. The net energy flux is upward. In the range 3-10 mHz, a full FOV space and time averaged flux of 2700 W m$^{-2}$ (lower layer 80-130 km) and 2000 W m$^{-2}$ (upper layer 130-180 km) is concentrated in less than 1% of the solar surface in the form of narrow (0.3$"$) AE. Their total duration (including rise and decay) is of the order of $10^{3}$ s. Inside each AE, the mean flux is $1.6 10^{5}$ W m$^{-2}$ (lower layer) and $1.2 10^{5}$ W m$^{-2}$ (upper). Each event carries an average energy (flux integrated over space and time) of $2.5 10^{19}$ J (lower layer) to $1.9 10^{19}$ J (upper). More than $10^{6}$ events could exist permanently on the Sun, with a birth and decay rate of 3500 s$^{-1}$. Most events occur in intergranular lanes, downward velocity regions, and areas of converging motions.Comment: 18 pages, 10 figure

Space-Time Distribution of G-Band and Ca II H-Line Intensity Oscillations in Hinode/SOT-FG Observations

We study the space-time distributions of intensity fluctuations in 2 - 3 hour sequences of multi-spectral, high-resolution, high-cadence broad-band filtergram images (BFI) made by the SOT-FG system aboard the Hinode spacecraft. In the frequency range 5.5 < f < 8.0 mHz both G-band and Ca II H-line oscillations are suppressed in the presence of magnetic fields, but the suppression disappears for f > 10 mHz. By looking at G-band frequencies above 10 mHz we find that the oscillatory power, both at these frequencies and at lower frequencies too, lies in a mesh pattern with cell scale 2 - 3 Mm, clearly larger than normal granulation, and with correlation times on the order of hours. The mesh pattern lies in the dark lanes between stable cells found in time-integrated G-band intensity images. It also underlies part of the bright pattern in time-integrated H-line emission. This discovery may reflect dynamical constraints on the sizes of rising granular convection cells together with the turbulence created in strong intercellular downflows.Comment: 24 pages, 15 figure

Solar Magnetic Carpet I: Simulation of Synthetic Magnetograms

This paper describes a new 2D model for the photospheric evolution of the magnetic carpet. It is the first in a series of papers working towards constructing a realistic 3D non-potential model for the interaction of small-scale solar magnetic fields. In the model, the basic evolution of the magnetic elements is governed by a supergranular flow profile. In addition, magnetic elements may evolve through the processes of emergence, cancellation, coalescence and fragmentation. Model parameters for the emergence of bipoles are based upon the results of observational studies. Using this model, several simulations are considered, where the range of flux with which bipoles may emerge is varied. In all cases the model quickly reaches a steady state where the rates of emergence and cancellation balance. Analysis of the resulting magnetic field shows that we reproduce observed quantities such as the flux distribution, mean field, cancellation rates, photospheric recycle time and a magnetic network. As expected, the simulation matches observations more closely when a larger, and consequently more realistic, range of emerging flux values is allowed (4e16 - 1e19 Mx). The model best reproduces the current observed properties of the magnetic carpet when we take the minimum absolute flux for emerging bipoles to be 4e16 Mx. In future, this 2D model will be used as an evolving photospheric boundary condition for 3D non-potential modeling.Comment: 33 pages, 16 figures, 5 gif movies included: movies may be viewed at http://www-solar.mcs.st-and.ac.uk/~karen/movies_paper1

Can we Determine Electric Fields and Poynting Fluxes from Vector Magnetograms and Doppler Measurements?

The availability of vector magnetogram sequences with sufficient accuracy and cadence to estimate the time derivative of the magnetic field allows us to use Faraday's law to find an approximate solution for the electric field in the photosphere, using a Poloidal-Toroidal Decomposition (PTD) of the magnetic field and its partial time derivative. Without additional information, however, the electric field found from this technique is under-determined -- Faraday's law provides no information about the electric field that can be derived the gradient of a scalar potential. Here, we show how additional information in the form of line-of-sight Doppler flow measurements, and motions transverse to the line-of-sight determined with ad-hoc methods such as local correlation tracking, can be combined with the PTD solutions to provide much more accurate solutions for the solar electric field, and therefore the Poynting flux of electromagnetic energy in the solar photosphere. Reliable, accurate maps of the Poynting flux are essential for quantitative studies of the buildup of magnetic energy before flares and coronal mass ejections.Comment: Solar Physics, in press. 14 pages, 3 figure

The Structure and Dynamics of the Upper Chromosphere and Lower Transition Region as Revealed by the Subarcsecond VAULT Observations

The Very high Angular resolution ULtraviolet Telescope (VAULT) is a sounding rocket payload built to study the crucial interface between the solar chromosphere and the corona by observing the strongest line in the solar spectrum, the Ly-a line at 1216 {\AA}. In two flights, VAULT succeeded in obtaining the first ever sub-arcsecond (0.5") images of this region with high sensitivity and cadence. Detailed analyses of those observations have contributed significantly to new ideas about the nature of the transition region. Here, we present a broad overview of the Ly-a atmosphere as revealed by the VAULT observations, and bring together past results and new analyses from the second VAULT flight to create a synthesis of our current knowledge of the high-resolution Ly-a Sun. We hope that this work will serve as a good reference for the design of upcoming Ly-a telescopes and observing plans.Comment: 28 pages, 11 figure

The Wave-Front Correction System for the Sunrise Balloon-Borne Solar Observatory

This paper describes the wave-front correction system developed for the Sunrise balloon telescope, and it provides information about its in-flight performance. For the correction of low-order aberrations, a Correlating Wave-Front Sensor (CWS) was used. It consisted of a six-element Shack - Hartmann wave-front sensor (WFS), a fast tip-tilt mirror for the compensation of image motion, and an active telescope secondary mirror for focus correction. The CWS delivered a stabilized image with a precision of 0.04 arcsec (rms), whenever the coarse pointing was better than ± 45 arcsec peak-to-peak. The automatic focus adjustment maintained a focus stability of 0.01 waves in the focal plane of the CWS. During the 5.5 day flight, good image quality and stability were achieved during 33 hours, containing 45 sequences, which lasted between 10 and 45 min. © 2010 The Author(s)

Modeling the Subsurface Structure of Sunspots

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.Comment: 73 pages, 19 figures, accepted by Solar Physic